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SUMMARY 

In many optimal control problems, it is reasonable to neglect certain 
of the terms in the system equations which are thought to have small effects 
in order to make the solution of the problem more tractable. When this 
results in a decrease in system order, however, the resulting approximation 
will not in general be able to satisfy all of. the system boundary conditions 
and thus the approximation will not be valid, at least in a small region. 

In this case, the system is said to be singularly perturbed. Unfortunately 
most of the results of singular perturbation theory to date have been 
concerned with initial value problems whereas optimal control problems are 
of two-point boundary value type. The portions of this theory applicable 
to the open loop state regulator problem are reviewed in this paper. 

For obtaining approximate solutions to the state regulator problem 
the method of matched asymptotic expansions is employed. This method has 
been developed in connection with certain fluid mechanics problems and 
is applicable to nonlinear as well as linear problems. It has been found 
in the past to be advantageous not to formulate this method generally but 
to apply it to each individual problem and this approach is adopted here. 

A general recipe for the method is given and its application is illustrated 
by using the method to obtain an approximate solution to a simple, specific 
state regulator problem; 
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I. INTRODUCTION 

In the state regulator problem, it is desired to keep the state of 
a controlled, linear, dynamic system near zero without using excessive 
control expenditure. It is well known that application of the necessary 
conditions for optimal control to the state regulator problem results in 
a two point boundary value problem for a system of linear ordinary differential 
equations. This paper is concerned with obtaining approximate solutions 
to such systems of equations by using singular perturbation methods. 

Although the feedback solution of the state regulator problem is known^ , 
we do not take advantage of this fact since we desire a method which is 
applicable to nonlinear as well as linear systems. 

Since efficient algorithms are available for the solution of systems 
of linear equations, the question arises as to why methods of obtaining 
approximate solutions need be studied. Approximate methods may be desirable 
or even necessary when computational speed or storage capacity is at a ' 
premium; for example in preliminary design in which a large number of trial 
systems must be evaluated or for use in a system on board a flight vehicle. 
Another important reason for studying linear systems is to gain insight 
into techniques which may be later applied to nonlinear systems. The primary 
factor determining computational requirements for solving state regulator 
problems is system dimensionality. Thus approximate solutions based on 
reduced order systems are often employed in practice. 

A standard technique of obtaining approximate solutions of mathematical 
problems is to introduce perturbations about a nominal solution. This 
technique is particularly useful in problems in which there is a "small 
parameter" present because in this case the nominal solution and the method 
of introducing the perturbations are suggested in an obvious way. When 
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the nominal solution is obtained from a reduced order system, it will not 

in general meet all of the boundary conditions of the complete system and 

thus the nominal solution will not approximate the solution of the complete 

system uniformly. This loss of uniform convergence characterizes singular 

perturbation theory. This theory is reviewed in references 2 and 3. 

The purpose of this paper is to briefly review the singular perturbation 

theory of ordinary differential equations as it applies to linear systems 

with two point boundary values and to illustrate application of this theory 

to optimal state regulator problems by solving a simple example. In this 

example, the theory is implemented by using the technique of matched 

4 5 

asymptotic expansions ’ which has been used to solve boundary layer 
problems in fluid mechanics. Use of singular perturbations and matched 

asymptotic expansions in optimal control problems was suggested by 

6 7 8 

Kelly ’ , and has been studied most extensively by Hadlock . 

In the recent literature, there have been two distinct methods (other 

than that under consideration here) proposed to obtain approximate solutions 

based on reduced order models of linear systems. In the first of these, 

the existence of a small parameter is assumed and the control is expanded 

« 9 10 

in a Taylor series in the small parameter and the first two terms retained ’ 
The zero order term is the optimal control of the reduced order system 
(obtained by setting the small parameter equal to zero) and the first 
order term is used to account for the initial conditions. In the second 
method, the characteristics of the model, including its order, are adjusted 
to give the best agreement (integral norm) between the output of the model 
and the output of the system using the model optimum control 11 . Either 
or both of these methods suffer from the following deficiencies: 

(1) inadequate treatment of control variables, (2) lack of proof of uniform 
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convergence to exact solution, (3) failure to deal directly with loss 
of boundary conditions, and (4) inability to generate a sequence of successively 
better approximations. The method of singular perturbations, using matched 
asymptotic expansions, does not suffer from these shortcomings. Further, 
as opposed to the methods of references 9, 10 and 11, it is applicable 
to nonlinear systems. 
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II. STATE REGULATOR PROBLEM 


In this section the solution of the conventional state regulator 
problem is first presented. The analysis follows that of reference 1, 
but with the nomenclature and sign conventions of reference 12, and is 
only summarized in this paper. Next, the problem is reformulated, in a 
form suitable for application of singular perturbation theory, for systems 
in which there is a small parameter present. 


Conventional Problem 

The most general problem to be considered is as follows. It is desired 
to minimize the cost functional 


T 

J = \ S [ < X (t)» £. (t) X (t) > + < a (t), R (t) a (t) >] dt (2.1) 


for the following system 


^ (t) = A (t) )L (t) + B_ (t) a_ (t) ; x. (0) = Xq (2.2) 


where x^ is an n dimensional state vector, a_ is an r dimensional control 

vector, and <•, •> is the scalar product on R n . It is assumed throughout 

that (a) T is fixed, (b) £ is positive semi -definite on 0*t*T, (c) is 

positive definite on 0^t<T, (d) optimal control exists, and (e) ' c_(t) is 

unconstrained. To formulate the necessary conditions for optimal control, 
12 

the 3C function is formed as* 


ft=-i<x, £x 


1 

2 < — * 


,B.£. > + < Ax_>A_ > + < j^a.>X_ )> 


(2.3) 


♦Hereafter functional dependence will be omitted when this does not result 
in confusion. 
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\ 


where A has been taken as -1. The adjoint vector A is defined as 
o — 

the solution of 

A = - 1? ; A (T) = 0 (2.4) 

- 3X - - 

Since QT = Q, this gives 

i = £ x - A" A_ (2.5) 

The optimal control ^ is that which maximizes since is positive 
definite ( R/ = jR, |&| +0) this gives 

gK i i 

si - - \ a £ - \ r° + s' a - a 

o = R- 1 B'i (2.6) 

g2 -i/ 

That (2.6) maximizes 3C follows from — - - R and the positive definiteness 

oO dO = 

of R. Putting (2.6) in (2.2) and adding (2.5) gives a 2n dimensional linear 
homogeneous system with evenly split boundary values 

x (t) = A (t) x. (t) + £ (t) X (t) ; x (0) = x^, 

A. (t) = £ (t) x (t) - A' (t) 1 (t) ; x (T) = 0 ( 2 . 7 ) 

where 

S = B R’ 1 B ' (2.8) 


*The adjoint vector used here is the negative of the costate vector used 
in reference 1; thus 3C is maximized to minimize J, whereas it is minimized 
in reference 1. A further difference from reference 1 is in the handling 
of the integrand of the cost functional. 
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Integration of this system gives the (open loop) solution of the problem. 

It can be established for the system (2.7) that there exists a symmetric 
matrix j£(t) such that 

± = K x (2.9) 

and this relation may be used to eliminate A in the first of (2.7): 

A = (4 + 1= K) A ^ A = Ao (2.10) 

The second of (2.7) is used to obtain a relation for J<; differentiating (2.9) 
and using (2.5) and (2.10), 

• • 

i = K x + K x 

£x-A^Kx=Kx+K(A+SK)x 
Since A 'fn*, this implies 

K = - A' K- KA - KU ; K (T) = 0 (2.11) 

where the boundary condition follows from (2.9) and \ (T) = £ provided 
_x (T) =f= (L Putting (2.9) in (2.6) gives the closed loop (feedback) control 
as 

a = R' 1 B' K x (2.12) 

Thus solution of the nonlindar (Riccati) system (2.11) and then integration 
of (2.10) gives the solution to the original problem; this is an alternate 
procedure to solving (2.7) and has the advantage of giving closed loop control. 


*The trivial case Xq = £ is excluded. 
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Problem With a Small Parameter 

In many systems, there is present a parameter which may be identified 

as "small" based on physical reasoning. For example, in viscous fluid 

flow, the viscous effects are small in many problems, at least in large 

regions of the flow field. Thus the viscosity is often treated as a small 

4 5 

parameter in the Navier-Stokes equations * . In other problems, however, 
no such parameter appears from physical grounds and a parameter may be 
artificially inserted to suppress terms in the equations which are expected 
to have relatively small effects. There is nothing new in this procedure, 
since this is in effect what the analyst does everytime he formulates a 
problem in dynamics, i.e., decides which effects to account for and which 
ones to ignore. In a flight mechanics problem for a manned vehicle, for 
example, a complete set of equations of motion would consist of the coupled 
system of the six equations of rigid body motion of the vehicle as a whole, 
the equations describing the dynamics of the control systems, the equations 
describing the dynamics of the pilot's arm and foot, etc. It is obvious 
that many of these effects can be neglected if, say, the vehicle trajectory 
is the only thing of interest. In what follows, it is assumed that the 
small parameter in the system has either been identified or inserted in 
a suitable manner. 

Consider the problem of minimizing 

1 C T 

J = 2 Jq C < 1 ( e, t) , Q ( e, t) x ( e , t) > + <a ( e» t) , R ( e, t) a ( e, t) >] dt 

(2.13) 


where the system is given by 


4 (e) X (e, t) = A (e, t) X (e, t) + B ( e , t) (e, t) ; X (0) = X Q (2.14) 


with the same assumptions as before, and where e is a scalar parameter. 
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Without loss of generality it may be assumed that ^ is a diagonal matrix, 
|£|fO> and that all matrices in (2.13) and (2.14) are of order zero or 
higher in the parameter e, that is, for example. 


lim = £ 

e+0 


We call this the singularly perturbed state regulator problem. The most 
common form of 4 in applications is a matrix with the first m elements 
equal to unity and the last n-m equal to e , i.e. 


e = 


1 


\ 


•1 


0 


e / 


(2.15) 


Multiplying (2.14) by £ results in 


x = £ 1 A x + e“ 1 B a 


- l 


(2.16) 


to which the results from the conventional system may be applied. Noting 
that, for example, (£ _1 A) ' = A ' jf 1 , (2.6), (2.7), (2.8) give 


x = £ _1 A _x + £ 1 S 4 " 1 A ; x_ (0) = 

\ = g x - A/ £ 1 1 ; X CO = 0 (2.17) 

wi th 

a = R “ 1 Bl ' 1 A (2.18) 

where the adjoint vector is denoted by A, Making the change of variable 

4 = 4 ^ (2.19) 

gi ves 

£ ( e) X (e, t) = A (e, t) X (e, t) + S (e, t) _A (e, t) ; X (0) = X^ 

£ (e) A. (e, t) = Q (e, t) X (e, t) - A/ (e, t) _A (e, t) *, A. (T) =0 


( 2 . 20 ) 
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with a_ and i given by (2.6) and (2.8). Equations (2.20) may be called 
a singularly perturbed linear system; integration of this system gives 
the open loop solution to the problem. 

The closed loop solution may be obtained in the same manner as for 
the conventional system. Multiplying (2.9) by jl, differentiating, and 
using (2.20) gives 

ii = Q - A' K - 4 U ' 1 A - iKi 1 S K ; K (e, T) = 0 ( 2 . 21 ) 

which may be termed a singularly perturbed Riccatti differential equation. 
After (2.21) has been solved for K> the state history may be found from 

& A = (A + II ) x ; x (0) = 5^ (2.22) 

and the feedback control from (2.12). An important special class of systems 
are those for which (a) the system is controllable, (b) A, B^, Q, R are 
time-invariant, and (c) T ->- «. In this case, K is not a function of time 1 
and (2.21) reduces to an algebraic equation 

0 = Q - A' | - £ | i 1 A - iKi 1 IK (2.23) 

The analyses of references 9, 10 and 11 utilize the closed loop 
solution to the problem as a starting point. In this paper, we use the 
system (2.20) as a starting point since we desire a general technique; 
closed loop solutions are not generally available for nonlinear systems. 
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III. SINGULAR PERTURBATION THEORY AND THE 
METHOD OF MATCHED ASYMPTOTIC EXPANSION 

It is desired to review results of singular perturbation theory as 
it applies to linear, two-point boundary value problems and to develop 
the method of matched asymptotic expansions for the approximate solution 
of such systems. Specifically, the system (2.20) is of interest; to 
formulate this system more concisely, let 



The resulting problem is called P £ : 

P E : E (e) W (e, t) = C (e , t) W (e, t) ; R (t) W (t, 0) + S ( e ) w (e, T)=k(e) 

(3.2) 

where the boundary conditions have been generalized by ..a 1.1 owing "mixing" 
and dependence on e. The reduced (or degenerate) problem associated with 

(3.2), P , is obtained by setting e = 0: 

P Q : E (0) w (0, t) = C (0, t) w (0, t) ; R (0) w (0, 0) + R (0) w (0, T) = k (0) 

(3.3) 

Of particular importance is the case when 4 has the form given by 


(2.15). Let 
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x = 




z 



A = 


hi kl 

; S = 

Si h. 

; Q = 

£l ^2 

hs hk 


5=3 hk 


0.3 0.4 


(3.4) 


(3.5) 


where and £ have dimension m, z. and £ have dimension n-m, matrices 
subscripted 1 are mxm, and those subscripted 4 are (n-m)x(n-m). Then the 
system (2.20) assumes the following form, called the P £ problem. 


£=A 1 x+A 2 z + Si£+S 2 a.;y_(0) = ^ 

4 = A 3 l + A^ z + S 3 p. + ; z (0) = ^ 

P £ : (3.6) 

A = 9,i H + O 2 1 “ A{ £ - A 3 £ ; £ (T) = 0 

4 = 9,3 1 + 0.4 l - hk £ - hk a ; i CO = o 


As before, the reduced problem associated with (3.6) is denoted P q : 

l=bj¥_+t2L + h£- + hZ>X-W = ^ 

0 = A 3 ^+A 4 z + S 3 £.+ Sj + g_ 

P o : (3.7) 

A = Qi £ + Q 2 L ~ hi R " 45 R » R (T) = £ 

o = Q 3 l + % L - hi R - ^4 a 

where the matrices are evaluated at e = 0. 

An even more concise statement of this system is obtained as follows; 


let 



- 13 - 



F = 


4i 

Qi -a: 






(3.8) 


(3.9) 


where and are 2m and 2n-2m dimensionaly , respectively; then (3.6) may 
be written: 

d jj (e , t) 

= £ (e, t) U_ (e, t) + D (e, t) V. (e, t) 

dt 

u i ( e> °) = x o * 1 = 1s m 
i 

iif (e , T) =0 ; i = m + 1 , " , 2m 

P e : (3.10) 

d V (e , t) 

e = jj ( e , t) U (e , t) + G ( e , t) V ( e , t) 

dt 


v. (e , 0) = x ; i = 1, ", n-m 
i+m 

v. (e , T) =0 ; i = n - m + 1, 2n-2m 

with associated reduced problem 
d (0, t) 

= £ (0, t) u (0, t) + D ( 0, t) v (0, t) 

dt 

0 = H (0, t) u (0, t) + G (0, t) v (0, t) 

P : 

Uf (0, 0) = x Q ; i = 1, ", m 
u. (0, T) =0 ; i = m + 1, ", 2m 


(3.11) 
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Relation of Solution of Full Problem to That of Reduced Problem 

If e is a small parameter, it is natural to attempt to approximate 

the solution to or to P e by solving the associated P Q *. However, a 

fundamental difficulty is encountered at once, since the P solution will 

o 

in general satisfy only 2m of the 2n boundary conditions of P £ or P £ , 
and thus convergence of the solutions of the reduced problems to the 
solutions of p£ and P £ will be non-uniform at best. 

Unfortunately, almost all of the theoretical work in singular perturba- 
tion theory has dealt with either boundary value problems for finite order 
scalar differential equations or with initial value problems for first 
order systems, either linear or nonlinear; this work is reviewed in refs. 

2 and 3. Two authors have, however, considered two point boundary value 

O 10 

problems for first order systems and the following result gives the 
relation of the solution of P„ to that of P . 

Theorem 1: For P £ , suppose that there exists an e 0 such that on the 
domain D = { (e , t) | 0 1 e le , 0 £ t IT} we have 1) A, j|, Q, jl are continuous 
in both arguments, 2) £ is nonsingular, and 3) n-m of the eigenvalues of 
£ have negative real parts. Then, 1) a unique solution to P £ exists, denoted 
iu (e, t), v. (s » t), 2) a unique solution to P Q exists, denoted £ (0, t), 

X (0, t), 3) lim u_ (e, t) = u_ (0, t) uniformly on 0 it IT, and 

£->0 

4) lim v (e, t) = v. (0, t) uniformly on any closed subinterval of 
0 <t <T. Hypotheses 1 and 2 are sufficient to arrive at conclusions 1 
and 2; 3 is a stability condition and its significance will become 
apparent subsequently. 

13 

The theorem may be generalized in several important ways. Harris 
has considered systems of the type P^- except that his form of ^ is somewhat 

*In what follows, the context will make it clear whether (3.3) or (3.11) 
are to be used for P . 




- 15 - 


more restricted than that in P £ . It is found that several complications 

occur. First, with boundary conditions of the mixed type as in it is 

not immediately clear which boundary conditions should be cancelled for 

3 13 

the reduced problem P . Thus a "cancellation law" ’ must be formulated 
and an additional hypothesis added to the theorem to the effect that this 
law is satisfied unambiguously. Also, hypothesis 2 of Theorem 1 must be 
reformulated in such a way that all matrices which must be inverted to 
solve P Q are nonsingular. Another important generalization is to the time 
independent case discussed previously for which T-* 00 . This case has 
been studied by Hoppensteadt^ for the initial value problem and a result 
similar to Theorem 1 has been obtained. 

For nonlinear problems, extensive modification of the theorem is 
required. In this case, it must be assumed that there exists an isolated 
root of the reduced system (3.11b), that this root is asymptotically stable 
in the sense of Liapunov, and that the boundary conditions are in the domain 
of influence of this root. For nonlinear initial value problems this theory 
was developed by Tihonov and is reviewed in ref. 2. An interesting geometrical 
interpretation of Tihonov's theory is contained in ref. 15 where it is seen 
to include the phenomenon of relaxation oscillation. A preliminary extension 
of these results to the two-point boundary value problem is contained in 
ref. 8. 

For the feedback solution to the state regulator problem, there is 
a result analogous to Theorem 1; under suitable hypotheses it can be shown 
that the limiting form of the gain matrix J< (e, t) obtained by setting e = 0 
in (2.21) is the same as the gain matrix of the reduced system. 
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Boundary Layer Equations 

Theorem 1 indicates that v^ ( e ,t) will change very rapidly from its 
boundary values at t = 0 and t = T given by (3.10) to values approximating 
v^ (0, 0) and v_ (0 , T), respectively. These rapid changes are called "boundary 
layers," a term which comes from fluid mechanics in connection with certain 
problems of viscous flow past a solid body. In these problems, the viscosity 
is a small parameter multiplying the highest derivatives in the Navier-Stokes 
equations. If this parameter is set to zero, the hydrodynamic system of 
equations results (reduced system); the solution of this system violates 
the no-slip boundary condition at the body surface. Thus in a thin layer 
of fluid near the surface of the body, the boundary layer, the velocity 
varies rapidly from zero on the surface of the body to the value given by 
the hydrodynamic solution. The boundary layer is characterized by high 
velocity gradients and high viscous forces. It is interesting to note that 
if the hydrodynamic solution is used without accounting for the boundary 
layer, the erroneous result that there is no streamwise force (D'Alembert's 
paradox) is obtained, even though the flow field is accurately described in 
all but a small region. This is indicative of the caution which must be 
exercised when drawing conclusions form reduced order systems. 

The phenomenon of boundary layers occurs in all singular perturbation 
problems. For the system (3.10) there will be boundary layers at both 
t = 0 and t = T; in these layers v^ ( e , t) will change rapidly from its "inviscid 
value _v (0, t) to its boundary values (e, 0) and v_ (e, T) in a time of order e 
To study this rapid change, we "stretch" the time scale by introducing the 
transformations 



- 17 - 


T _ _ 


t T 


T - t 


Using (3.12) in (3.6) the boundary layer equations associated with P 
at the initial point are found to be 

% 'b 'b 'b 

h = e A x & + E A2 Z 2 + e Sj £ 2 + e £ 2 ; & ( 0 ) = ^ 

'V % 'V 

li = A 3 Xi + A, Zj + s 3 + S 4 ; z 2 (0) = 2 ^ 

Pni • 'b 'b 'Xj 'b ^ 

1 £.1 = e £i y.1 + s £2 L\ ~ e A 2 - e A 3 




'b 


(3.12) 


(3.13) 


ii = £3 Zi + ^4 1 1 - £2 £1 - g . 2 

d Zi 

In these equations, the independent variable is T , for example y_ x = ^ — 

'b 0 T 0 

and Aj ( e , Tq ) = A 2 ( e , ex 0 )- The boundary layer system at the final point, 
Pgi^s is same as (3.13) except that the boundary conditions are replaced 
by £ 2 ( 0 ) = 0^, £2 (0) = 5. and The independent variable is T j. The abbreviated 
forms of these equations are 


d £.1 (e 5 x„) 


dx. 


e Fj (e> x q ) Uj (e> x q ) + e £j (e, x q ) Xl (e>T 0 ) 


U-, (e, 0) = x ; i = 1, ", m 

i i 


BL 1 ■ 


d Zl (e >x_) 


(3.14) 


dx. 


Hi (e» X ) Uj (e,x ) + G 2 ( e , Xq ) Vj (e,x ) 


V, (e, 0 ) = X 
i °i+m 


i = 1, n-m 
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d U2 ( £ » Tj) 
dx T 


£ £2 ( £> T y) M2 ( £ » T j) + e £2 ( £s T y) M2 (e » Ty ) 



u 2 ( e . 0) = 0 ; i = m + 1 , •••, 2m 
i 

d M 2 (e , tj) 

= if? ( £ > T j) M2 ( £ > T j) + J2 ( £ » t t) M2 ( £ > Tt ) 

dtj 1 1 1 


(3.15) 


v 2 (e , 0) =0 ; i = n - m + 1 , - , 2n - 2m 
i 

where, for example, fi (e , t ) = £ (e , eT Q ) and ( e , t t ) = F ( e , T - ex-p) . 


Matched Asymptotic Expansions 

It is natural to attempt the solution of (3.10) by asymptotic expansion 

in the small parameter e. Asymptotic methods for both regular and singular 

problems are reviewed in refs. 4 and 17. Before proceeding, some definitions 

are needed. Given functions f (e ) and g(e), we write f( E ) = 0 (g(e)) if f(e)/ 

Ml 0 as e -*■ 0 + . 
e; 

The sequence f n (e) is an asymptotic sequence if f^ (e) = 0 (f^ (e)) as 

e ■> 0 + ; f- (e) is an asymptotic expansion of the function F(e) if it is 

an asymptotic sequence and if 
n 

I F < e > - ft) Vl< £ > I ■ 0 (f n < E » 

which implies 
n 

l p U) - | 0 » i f i (0 I =°(f ntl (c)) 

We call ^ a^f^e) the nth order approximation of F(e). It should be 
i=0 

noted that different asymptotic sequences may lead to different asymptotic 
expansions for the same function but that the expansion in terms of a given 


g(e) is bounded as e-*0 + and f(e) = 0 (g(e)) if 
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sequence is unique. If the sequence is infinite, we write 

CO 

F(e) 'v, X a.f.(e) 
i=0 

and say that the series a^f.fe) is asymptotical ly convergent to F(e). The 

series need not be convergent in the ordinary sense. The only asymptotic 

sequence which will be used in this paper is the power series, e 1 . 

There is at present no theoretically justifiable technique of obtaining 

an asymptotic solution to (3.10). The most general result available is due 
1 8 

to Vasileva and deals with the Cauchy problem for nonlinear first order 
systems. In her method a solution is sought in the form of a sum of two 
asymptotic expansions, one which approximates the solution "inside" the 
boundary layer and one "outside," and an explicit method for obtaining the 
terms in these expansions is obtained. This method is not used in the present 
paper, but the following result of ref. 18 is of interest in motivating the 
use of asymptotic methods in singular perturbation problems. 

Theorem 2: Consider the system 

H Y 

= f ( x >y» e »t) ; x( o) = X Q 

e ^t = 9(x,y,e,t) ; y ( 0) = y Q 

with associated reduced system, 

= f(x,y,0,t) ; x(0) = x Q 


0 = g(x,y,o,t) 
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and associated boundary layer system (t = ex). 


dx 

dt 

ef (x,y,e,ex) 

• 

X 

o 

II 

X 

& = 
dt 

g(x,y,e,ex) 

*< 

o 

II 

<< 


and suppose that: 1) f and g are n + 1 times continuously differentiable 

in x, y, and e in some domain n and continuous in t for 0 < t s T, 2) the 
full system and the reduced system both have unique solutions in q, 3) there 
exists in isolated root y = <p (x,t) of the reduced problem 0 = g(x,y,0,t) in 
fi, and 4) the root y = tj> (x,t) is an asymptotically stable equilibrium point 
of the boundary layer system. Let the asymptotic sequences x^ (t)e 1 , "x - ( t ) e 1 , 
y^ (t)e 1 , y • ( t ) e 1 be those discussed in ref. 18. Then these sequences are 
asymptotic expansions of the solution the full system, x(e,t) and y(e,t), that 
is 

n n 

|x(e,t) - Z x,(t)s'' - 2 X.(t) E 1 | = o( E n ) 
i=l 1 i=l 1 

n n 

|y(e.t) - 2 y.- (t)e 1 - 2 = o(e n ) 

i=l 1 i=l 1 

in n for 0 x t T. 

Q 

Hadlock formally applies Vasileva's method to two-point boundary value 
problems and it appears that the method is applicable to this class of 
problems as well. The method to be used in this paper was originally developed 

in fluid mechanics and is called the method of matched asymptotic expansions 

(or the method of inner and outer expansions) ’ . The method was first conceived 
by Prandtl and has been formalized by Kaplun, Lagerstrom and others. In this 
method it has been found to be advantageous not to give an explicit general 
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formulation of the required expansions but to formulate only a general recipe 
and then treat each new problem individually, and this approach will be used 
here. This is due to the large diversity of behavior which may be encountered 
in even the most elementary problems as well to the great algebraic complexity 
generally involved. Theoretical justification of the method is not of great 
importance since it becomes readily apparent in the course of a solution whether 
or not the method is working. However, it seems likely that t lie hypotheses of 
Theorem 1 plus increased restrictions on the smoothness of the system matrices 
in the domain of interest are required. Specifically, hypothesis 1) will be 
replaced by: A, B = , £, Rare of class C°° in t and have asymptotic expansions in 
e in domain D. 

In the method of matched asymptotic expansions we proceed as follows: 

1) solve the "outer" system (3.10) asymptotically, leaving the boundary 
conditions free, 2) solve the "inner" or boundary layer systems (3.14) and 
(3.15) asymptotically subject to the appropriate boundary conditions, and 3) 
determine the remaining constants of integration by invoking a "matching 
principl e. " 

To solve the outer system set 


u[( e ,t) <v ^ i/ 1 ^ (t)e 1 
i=0 

CO 

F(e,t) -v 2 F (l) (t) £l 
i=0 

oo 

O. 2 H (1 W 
i=0 


co 



v 

to , 

I (t)c 

; v(e,t) 

^ Z 


i=0 

oo 


; D(e,t) 




i=0 


; G(e,t) 




i=0 



V 


(3.16) 


H(e,t) 
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Putting (3.16) in (3.10) gives 



i=0 dt 


i 


e 



CO 



+ 






1 


i=0 



dt 


i+1 

e 


CO 00 oo 



Equating the coefficients of successive orders of e of these series in the 
usual way results in a sequence of linear problems, each of which is of 
order 2m and hence 2m arbitrary constants will appear at each order. If it is 
assumed that is nonsingular these problems will have unique solutions. Note 

that the zeroth order problem is homogeneous while the higher order ones are 
inhomogeneous. 

To solve P DI , set 
BL 1 


^1 (e » T o> 





(3.18) 


in (3.14) to get 



(3.19) 
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where, for example. 



The zeroth order terms in these expansions must satisfy the initial conditions 
since these conditions are of zeroth order in e and the higher order terms 
must vanish at t = 0, that is 



x„ if i = 0 

°J 

0 if 1 > 1 


j = 1 


m 


v 


( 1 ) 



if i = 0 
if i > 1 


j = 1 


n-m 


(3.20) 


The system (3.19) is of ordinary perturbation type and has 2n constants of 
integration at each order; imposing the n conditions (3.20) leaves n arbitrary 
constants. Note that 2m of the equations (3.19) are simple quadratures. 
Similarly for Pg^ we set 


i!o(e»T T ) I ; v ? (e,t t ) ^ 2 v (^(t-jOe 1 (3.21) 

* 1 i=0 c 1 c 1 . i=0 c 1 


so that 


co 


2 




i 


€ 


% 



oo 


CO 


+ 
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subject to 

u 2 ^(0) = 0 ; i >_ 0, j = m+1 , . . , 2m 

j 

v 2 (i) (0) = 0 ; i _> 0, j = n-m+1 , 2n-2m (3.23) 

j 

where, for example, 

co 

F,U.t t ) % J F.WlOe' 

1=0 c 1 


The unknown constants in the outer and boundary layer expansions are 

determined by a matching principle. This principle assumes that there 

exist two overlap regions, one near t=0 and one near t=T such that in the 

first of these both the outer and the expansions are valid and agree, 

and in the second of these both the outer and P DI expansions are valid and 

[5 L ^ 

agree. In the limit, this leads to the matching conditions: 


U? u (i) (t) 


,im u . (1 >(t ) 

T -K» — 1 ' 0 

0 

,im V (’•)(, ) 

11m U (1) (T ) 
T ^2 l T J 


lim V (i) (T ) 

T —2 ' T ; 


(3.24) 
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These conditions are to hold for all orders i of e. It is evident that 
asymptotic stability of the boundary layer solutions is required. If it is 
assumed that n-m of the eigenvalues of are negative this stability is 
insured since n-m of the n constants appearing at each order of the boundary 
layer solutions may be chosen to suppress the unstable solutions, leaving m 
undetermined constants. The m constants in P m , the m in P m and the 2m 

D L i D L 2 

of P £ are determined by the matching conditions. It will be shown later 
that the matching conditions associated with z and c[ are redundant. 

The representation of the solution obtained by the method outlined above 
has the undesirable feature that different representations are required for 
different values of the independent variable. It is desirable to obtain a 
single expansion valid uniformly for 0 < t < T and this may be accomplished 
by adding the two boundary layer expansions to the outer expansion and sub- 
tracting out the common parts, i.e., the terms which cancel out in the matching. 
This expansion is 



(3.25) 
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To summarize the results of this section, we state the following: 

Theorem 3: Consider the system P £ as given by (3.6) or (3.10) and suppose 

that there exists an e 0 such that on D = {(e,t) |0<_e<_e Q , 0 <t<J} we have 1) A, 

£, £ are of class C°° in t and have asymptotic expansions in e, 2) is non- 

singular, and 3) n-m of the eigenvalues of are negative. Then the method 
of matched asymptotic expansions leads to a unique asymptotic representation of 

the solution of P . 

e 

It should be noted that the theorem does not state that the resulting 
asymptotic sequences converge to the solution of P . The theorem will be proved 
explicitly for the zeroth order in the following section (extension to higher 
orders is similar) where it is also shown that the zeroth order outer problem 
is precisely the reduced problem P Q . 


Zeroth Order Approximation 

To illustrate the method of matched asymptotic expansions, the zeroth 
order approximation will be worked out in detail. Retention of only zeroth 
order terms in the expansions gives a sufficiently good approximation for 
many applications. From (3.17) 


du 


(o) 


F<°> u (o) ♦ D<°> V<°> 
dt = — = — 


0 


H ( 0 ) u ( 0 ) + 3 ( 0 ) v (o) 


(3.26) 
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The solution of (3.26) is, referring to (3.8), 
= u (o) = e (£ (o) - D G (o) 1 £ (o) )t c 


r,(o)i 


r v (°)i 

Z 

= v (0) = _ G (o)-l H (o) 


(o) 

_£ J 

£ (o) 



where £ is a vector of 2m constants of integration. 


From (3.19) the zeroth order term in the P D1 expansion satisfies 

BL 1 


du. (0) 

-fc = 0 ; u -j l °'(0) = x Q ; i=l, .., m 

o i i 


d\/ 

~ar 


( 0 ) 


1— =H (o) u < 0) +G (o) v, (o) ; v, <o) (0) - x ; 1=1, n-n 

i °i+m 


The solution to this system may be written 

(or 


*1 

V o) 

,£] (0) J 


±1 





^o)" 1 H (o) f 4 



where is the modal matrix of and 


h - r ^i T ° 

VS e 


(3.27) 


(3.28) 


(3.29) 


where the C-j are the n-m constants of integration associated with the 
eigenvalues of which have negative real parts. The solutions 
associated with eigenvalues which have positive or zero real parts have 
been suppressed by setting their associated constants of integration to zero. 
The constants C _ -j and are related by 

— o 

a/ 0 , (o) 


Og 


( 0 ) J . g(o)" 1 u (o). 




(3.30) 


so that only m of their components are independent. The zeroth order term 
in P gL satisfies 


d -?^ (n) 

= o ; u 2 ioj (0) = 0 ; i = m+1 , . . , 2m 


dv 


(o) 


-j— — = H (o) u 2 (°) + £ (o) v 2 ^ ; v 2 (o) (0) = 0 ; i = n-m+1 , 


(3.31) 

2n-2m 


with solution 



(3.32) 


where 
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with 

~z 2 <°>( oj 
0 _ 

First consider the matching of ^ 0 ^ at the initial point. From (3.24) 
lim (o) _ lim (o) 

t-K) *- T -x» JCl 

0 

z io) (0) ■ 

Similarly, the matching of £1°^ at the final condition gives 


h g (0) c 2 



(3.33) 


£ (o) (T) = 0 


These conditions together with (3.26) imply that the zeroth order outer 

solution is just the solution to the reduced problem (3.7). The constants 

C have now been determined. The matching of at t = 0 and at 

^ * 

t = T gives the constants £-| and C^: 



The matching conditions for z ^ and at the initial and final points are 
then satisfied identically as may be seen from (3.27b), (3.29b), and (3.30b). 
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An Alternate Procedure 

In the preceding development the necessary conditions were applied 
to the full system and the resulting two point boundary value problem was 
solved by the method of matched asymptotic expansions. It is natural to 
ask whether or not the reverse procedure is valid, that is, is it permissible 
to first expand the state equations asymptotically and then to apply the 
necessary conditions to each order of the expansions? This latter procedure 
is attractive since it may be expected to result in less algebraic manipulation 
for a given problem. To pursue this question, we partition (2.14) in the 
manner of (3.6) 

y_ = Ai y + A? z + B _ 1 a 

ei = A 3 Z + A 4 I + Jh o 

where ^ is m x r, £2 (n-m) x r and the other quantities are as defined 

by (3.4) and (3.5). The following result gives the conditions for equivalence 
of procedures for the zeroth order. 

Theorem 4: Let the first two hypotheses of Theorem 1 apply and further 

suppose that the matrix A 4 has maximum rank (i.e., rank n-m) in a 
neighborhood of the extremal. Then the reduced problem obtained by setting 
e = 0 in the necessary conditions of the full problem is the same as the 
problem obtained by setting e = 0 in the state equations and applying the 
necessary conditions to the result. 

The significance of the additional hypothesis is clear; in either 
procedure, if it is satisfied n-m of the components of ( z , g) may be eliminated 
from the problem. When e is set to zero in the state equations, _z becomes 
effectively a control variable and the two control variables z and g_ are constrained 
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by a set of state dependent algebraic equations. The adjoint variables as- 
sociated with the state variables whose derivatives are neglected correspond 
formally to the ordinary Lagrange multipliers arising in the maximization of 
the 'K function in the second procedure. The proof of the theorem follows 
from section 3.6 of ref. 12. A more restrictive equivalence of procedures 
result than Theorem 4 is given in ref. 8 where it is shown that a sufficient 
condition for equivalence is the non-singularity of Al+< If one takes the 
viewpoint that a control variable is a variaDle which may be changed instan- 
taneously*, it may be concluded that, for problems in which a small parameter 
is to be inserted, the parameter should be inserted in such a way as to make 
the relatively "fast" variables behave as control variables in low order 
approximations. 


*I.e., controls may be piecewise continuous or even just measurable. 



- 32 - 


IV. AN EXAMPLE 


Problem Formulation 

Consider a forced two degree of freedom spring-mass-damper constant 
coefficient system with specified initial conditions. 


ay + by + cy = a ; y(0) = y Q , y(0) = y Q 


It is desired to select the control a(t) to minimize 


1 00 

J = 7 f (qny 2 + q 2 2y 2 + 2q 12 yy + Ra 2 )dt 

J o 

/x x \ 

In state variable form (let x x = y , x 2 = y, x = ( x 1 ) 


Xi = Xc 


; (0) = x ! 0 


axe 


J 4 


bx 2 - cx l + a ; x 2 (o) = x 20 


oc 

l 


< X, Qx > + R a 2 


dt 


If b or c is "small" this is a regular perturbation problem and if a is 1 
it is a singular perturbation problem. Suppose a is small and q n = q 12 
Without loss of generality, let q 22 = q 2 , R = 1 and a = e so that 


Xi = x 2 


; x x (0) = x ! o 


ex 2 = - bx 2 - cx x + a 


J 



(q 2 x 2 2 + a 2 )dt 


smal 1 11 
= 0. 


> x 2 (o) = x 2 0 
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The system matrices are 



(° ") 

\0 q 2 / 


&■(’) 

From (2.8) 



From (3.5) and (3.9) 

4 i = ( 0 ) > A 2 = ( 1 ) ; A 3 = (-c) ; A 4 = ( - b ) 

Si = ( 0 ) ; £2 = ( 0 ) ; S 3 = ( 0 ) ; 5* = 0) 

0,1 = (0) ; O2 = (0) j £3 = (0) '» O4 = (^2) 



Since £ is nonsingular with eigenvalues + sj q 2 2 +b 2 and - J q 2 2 +b 2 , Theorems 

1 and 3 hold (G = G^), and it may be expected that the method of matched 
asymptotic expansions will work on this problem. Further, [A 4 , B^] = (-b 1) has maxi 
mum rank (in fact A 4 is nonsingular) so that Theorem 4 holds and the alternate 
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procedure may be used, that is we may set e = 0 in the state equations, treat 
x 2 and a as control variables, and apply necessary conditions to obtain the 
reduced problem. This equivalence may be verified by direct computation for 
this simple example. 

In the following, both the exact and zeroth order asymptotic solutions 
for the case of no damping (b = 0) will be obtained. Consideration of this 
case exhibits the essential elements of the method while minimizing algebraic 
complexity. Note that now 
A„= ( 0 ) 

a-(° ’) 

\q 2 0/ 

and [Ai* Jbl = (0 1) so that Theorems 1, 3 and 4 still hold. However^ the 
alternate procedure will not be used in the following. 


Exact Solution 

The system to be solved is (see 2.20) 


*i 

eX 2 

Xi 

eX 2 


Set 


x 2 

-Cxj + a 2 
C \ 2 

q 2 * 2 ■ 

.st 


; X! ( 0 ) 

x 2 ( 0 ) 
Xi(~) 
X 2 (“) 


Xio 

X 2 0 

0 

0 


Xi = Aie^’ 1 ', x 2 = A 2 e^, , x 2 _ A 4 e 


St 


.St 


to get 

A 2 S = A 2 
eA 2 S = -cAi + A 4 


A3S = cA 4 


eA 4 $ - q 2 A 2 - A 3 



- 35 - 


Let 


P2 - A^/A^ - S 

u 3 = ^3/^1 = § ( £ S 2 + c ) 


P4 - A^/A^ = eS + c 
The characteristic equation is 

e 2 s 4 + (2eC - q 2 )S 2 + c 2 = 0 


with solution 
S 


1,2, 3, 4 - 


+ v^--- + V— - — 

-V 0 IV . d Q 


o 2 

£e e 


4 g 4 e 3 


The exact solution is represented as 

«1 - A^ 1 ^ 1 * + A^'e 52 * ♦ a/”.* 3 * + a/ 4 ^ 4 * 


*2 1^2 


'DAjt^e 51 * + ^’a/ 2 ^ 2 * .♦ ♦ „ 2 < 4 \( 4 >e S4t 


Xj ■ uj^Aj 11 ^ 1 * ♦ „ 3 < 2 >A 1 < 2 >e S2t ♦ ^V^ 3 * ♦ u^A^e^ 


^2 ~ ^4 


V\Me Sit ♦ V4 < 2 >A 1 < 2 >e Szt ♦ ^V 3 ^ 3 * ♦ „ W*/ 4 ^ 4 * 


The A^^ are determined by the boundary conditions; this leads to 


( 1 ) 


a/ 2 ) = 0 


( 3 ) . 


(4) 


x 20 " u 2 x 10 
y (3) (4) 

v*2 ^2 
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A. 


( 4 ) _ 


'20 


~T3) 

M2 " ^2 


u (3) X 
u 2 * 10 

wr 


where 



For future purposes, we expand the solution for small e and retain 
only the first terms. For S^: 
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v/cu 


(4) 


\f^2 


■ W «- C ^ 


(4) 


Thus 


Al<3) “'^: (jb Xl ° + X20 


A< 4 >»x 


2 \ H 2 
10 


so that the zeroth order approximation is 

c 


\/<b 


x 10 e 


e 

r- W 

A 1 “ " c v /^2 X 10 e 




A 


2 





+ cx 


10 
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Terms with factor e 


Tq^/e t 


will be important only for small t (t = 0(e) or 


less) and appear only in the expressions for the variables (x£, X£) which lose 
their boundary conditions at e = 0. It may thus be anticipated that these 

terms are the zeroth order terms in the boundary layer expansions; and that 

- t 

terms with factor e are the zeroth order terms of the outer expansion. 

This will now be shown by obtaining these expansions. 


Solution by Matched Asymptotic Expansions 

The zeroth order outer system is just the reduced problem, as has been 
shown earlier. 


X 1 x 2 
X 0 0 


; x lo (°) = x 1Q 


A, = CA 9 

L o 


; A, (“) = 0 
x 0 


0 = - CX, + A, 
L 0 ^0 


0 = q 9 Xp - A, 


Solution of this simple system is 


X 1 Q *10 e 


2 0 ^2 10 


c v/^2 

t== e 
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A l„ " " C 'fil X 10 e 


J*2 


A 2 q °0 ' X 10 6 




These functions do not satisfy the initial conditions but do satisfy the 
terminal conditions and hence a boundary layer will be needed at the initial 
point only. 

The zeroth order boundary layer (inner) system is obtained by using 
transformation (3.12) and then setting e = 0: 


I 


= 0 


I 


= 0 


; x^(0) x lQ 


Xo ~ ~ C X, + Xn , Xp (0) Xp„ 

C 1 A I C 1 I u 


Ap dp Xp - A-. 

I I I .... 

where ( )' = ^ ^ . Of the four constants of integration for this system, two 
T o 

are determined by the boundary conditions, one is needed to suppress the 
unstable root of the characteristic equation, and the fourth, say K, is to be 
determined by the matching conditions. The solution is 


1 



X 


10 



T 

0 


"*■ (X£Q “ K) 
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A lj q 2 
X 2j = °I 

The inner and outer solutions are now matched according to (3.24): 

X 1 : x 10 = x 10 

X 2 : " x 10 = X 20 ' K 

A 1 • - c \ f ^2 X 10 = q 2 ^ x 20 ~ ^ 

X 2 * C x 10 = C X 10 

These conditions are satisfied unambiguously by the unique value 
K = x 20 + ^ x 10 

(In the previous section it was shown that the matching of and \ 2 is 
redundant.) Thus the zeroth order inner expansion is 


(x 2Q - K) 

= c x 1Q - e ^ 0 


x lj " x 10 


^2 


X i n + 


2 i 10 


yq^ x 10 + X 20 ) e 
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v 


k lj " ' C ^ 2 X 10 


v/o? 


A 0 a x c x in - ( c Xjq + \J C ) 2 ^20 ^ ® 


l 2j ' 


'10 


Finally, the zeroth order of the uniform expansion obtained according 
to (3.25) is 

c 


X l u X 10 e 




— t 

_ / c \ ~ t cx io y^2 

% = (/> 2 Xl0 + X20 ) 6 ~^2 6 


A l .. _ c \/^2 X 10 e 


n /^2 


v/^2 


A 2 m a U = " ( C x 10 + V^2 x 20 ) e 


« + c x e ^ 
X 10 e 


which agrees with the first term expansion of the exact solution. A useful 
approximation to the solution has thus been easily obtained. 

This example shows that application of the method of matched asymptotic 
expansions in effect "splits up" a given problem into relatively easier 
problems, and it is this characteristic which makes the method attractive for 
complex problems. The method may be continued to obtain first order and 
higher approximations. 


Numerical Example 

To gain insight into the ability of the uniform zeroth approximation 
of the example to approximate the exact solution, consider the specific case 
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c = i> i » x io 


X 20’ 


e - 


1 

10 


so that e is only one order of magnitude smaller than the other system 
parameters. For this case, the exact solution is 


- .2746 e~ 8 ' 873t + 1.2746 e~ 1A27t 


'10 


'10 


2.436 e ” 8 ' 873t - 1.436 e -1 * 127t 


'10 


o aic -8.873t , , -1 . 1 27 1 

- 2.436 e + 1.436 e 


the outer solution is 




the inner solution is 



1 
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2 e 


■lOt 


- 1 



2 e 


-lOt 


and the uniform solution is 




-e _t + 


2e 


-lOt 




2e 


-lOt 


The various solutions are plotted on Figures 1, 2, 3. Semi-log paper 
is used to amplify the behavior at small t (boundary layer). The figures show 
that agreement between the exact and the approximate 'solutions is best at 
very small t and at very large t and worst in the overlap region e <_ t <_ 10e. 

It appears that the inner solution is better than the uniform solution at small 
t and that the outer solution is better at large t. It is seen, however, that 
the uniform solution gives a reasonable approximation to the exact solution; 
for smaller orders of e this approximation would, of course, be improved. 

The cost is computed from 


J 



*2 + a )dt 
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Thus the approximate solutions are unconservative (underestimate the cost) 
with the outer solution giving a better value than the uniform one. 

A Special Case 

To illustrate that for singular perturbation problems a slight 
modification of the system may lead to dramatic differences in the solution 
characteristics, the special case of x^g = 0, c = 0 is now considered. The 
zeroth, first, and second order terms of the uniform asymptotic expansions will 
be obtained. In this case a^ = 0 so that the system to be solved is 

x 1 = x 2 ; Xj(o) a 0 

£X 2 — O , X 2 (0) — X 2 Q 

ea = q 2 x 2 ; a(«>) = 0 

The exact solution is 
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x 2 = x 20 e 


& 


- e 


To solve the outer problem asymptotically, to second order, set 


ct - (t)c 1 ; Xn ^ ^ 6^ ( t ) e 1 ; a ^ 2 a-(t)e^ 

1 i=0 1 * i=0 1 i=0 1 


in the system equations and retain only second order terms and lower: 


a o + a l E + a 2 e = 6 o + B l e + B 2 e 


• • 2 2 

B Q e + = o Q + a^e + o 


• • 2 2 

°o £ + a l e = q 2 B o + q 2 B l e + q 2 B 2 e 


To zeroth order: 


a o = B o 

0 = °o 
0 = q,0 


2 o 


a = C 
0 0 

y =» Bq = 0 
a o = 0 


where C is an undetermined constant. The first order terms are 
o 


a l = 8 l 
B o = a l 
°o = q 2 6 l 


“i = c i 

e 1 = 0 
= 0 
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The second order terms are 


a 2 = e 2 
= ct 2 

°1 ~ ^2^2 


o t ^ ^2 

s 2 = 0 
©2 = 0 


The outer solution to second order is then 


x 1 - C + C-i e + C«e 
i 0 o 1 2 


Xo 0 ; a = 0 


The boundary layer equations are 


Xj = ex 2 ; x^O) = 0 


Xo - o j Xo (o ) x 


20 


a ' = q 9 x 


2 2 


where ( )' = - 4 ^ and t = ex. To solve these to second order set 
dx 


OO CO CO 

-I zl 6 . ( x ) e ^ ; Xo X, 2 n^CxJe 1 ; a 'v, zL v^(x)e^ 

1 i=0 1 L i=0 1 i=0 1 


in the boundary layer equations and retain terms up to and including second 


order: 
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5 o + + S 2 e ^ = n o e + n i e ^ 

x 2 . , 2 
1q + n£e + n 2 e = v 0 + v l e + v 2 e 

V o' + V 1 E + v 2 tZ s «2 n o + W + V2 e2 


To zeroth order: 


V° 


v - = q.n 
o H 2 o 

« o (0) = 0 

V°> = X 20 


6 o = 0 


n = —r 
0 v/^2 


% v^2 T " v^2 

( K 0 + V^2 X 20^ e ' K o e 


\ft\o T “ v/Pp 

v o = < K o + v^2 x 20^ e + K o e 


The first order terms are 


6 i 


n l 
v „ 


^(o) 

\(o) 


"0 

V 1 

Vl 

0 

0 


1 

6 1 ~ q 2 


yqy x 

( K 0 + v /^2 x 20 ) ( e -1 ) 


- x/q? T 

+ K 0 (e v 2 -1) 


k i / v/^2 x - >/V 

"i = 7 q 7 \° _ e ' 

V 1' K 1 (e^^e"^ 


The second order terms are 


5 2 


Ho “ 


2 

& 2 (°) 

n 2 (°) 


n l 

v 2 

q 2 n 2 

0 

0 


. _ K i / A2 t ■ y^2 x " 
2 ' ^7 V / 

^2 / y^ 2 x - y q 2 t 

f] O r=P" I. G “ G 


y^2 v 


V 2 = K 2 \ e 


V^2 T + - v^2 
+ e 


5. 

"2 
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The boundary layer solution to second order is 


X 1 ~ ^ 1 £ + ^ 2 £ 


x 2j = n o + n l £ + n 2 £ 


°I = V o + V + V 


The constants C Q , C 2 , K qJ Kj, K 2 are now determined from matching. Since 
all orders of the boundary layer expansions must remain finite as t-x», 


K 0 ~ \J^2 x 20 * ^1 0 > ^ 2 0 


so that 


'-M ( 1 - e' ^ 


x 20 e 


-\ fi ~2 


= - To? x 




2 20 


Matching with the outer solution gives 

x 20 

C o ” 0 * C 1 = yqj ’ C 2 ” 0 


so that 


X = M c 

*0 ^2 


x ? = 0 
0 


°o = 0 
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The uniformly valid expansions are 

_ >^2 

1 - e e 


l l 


U 


20 
A 2 


X 2 U x 20 e 


v/q ” 2 


a,, - 


v/q? x on e 


U vh 2 ^20 


Thus the boundary layer expansions are the same as the uniform expansions 
and the outer expansion has contributed nothing. Also note that the second 
order terms are all zero so that higher order terms will be zero also, 
that is 


a. = 6- = a. = <5. = n. = v. =0 M i>2 
iiiiii - 

This is verified by noting that the uniform expansions above are precisely 
the exact solutions. It may be concluded that for this special case, the 
method of matched asymptotic expansions has not resulted in, computational 
simplification. 
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CONCLUDING REMARKS 

Singular perturbation theory and the method of matched asymptotic 
expansions have been applied to the state regulator problem of optimal control 
theory. Conditions have been stated under which the method is applicable. 

The method is illustrated by using it to obtain an approximate solution of a 
simple, specific, singularly perturbed state regulator problem. It is found 
that the method is easy to apply and results in a uniform approximation 
which gives good agreement with the exact solution in the entire range of 


interest. 
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Figure 2. Comparison of solutions for second state variable. 
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Figure 3. Comparison of solutions for control variable, 




